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Abstract 

The odd-even mass staggering in nuclei is analyzed in the context 
of self-consistent mean-field calculations, for spherical as well as for 
deformed nuclei. For these nuclei, the respective merits of the en- 
ergy differences and A^ 5 ) to extract both the pairing gap and the 



1 



time-reversal symmetry breaking effect at the same time are exten- 
sively discussed. The usual mass formula A®, is shown to contain 
additional mean-field contributions when realistic pairing is used in the 
calculation. A simple tool is proposed in order to remove time-reversal 
symmetry breaking effects from A^ 5 \ Extended comparisons with the 
odd-even mass staggering obtained in the zero pairing limit (schematic 
model and self-consistent calculations) show the non-perturbative con- 
tribution of pairing correlations on this observable. 
PACS: 21.10Dr; 21.10.Hw; 21.30.-x 

Keywords: Mean- field theories; Pairing correlations; Odd-even mass 
staggering; 
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1 Introduction 

The Odd-Even Staggering (OES) of binding energies is a common phe- 
nomenon of several finite many-fermion systems. In nuclei, it has been 
attributed to an experimental evidence of pairing correlations 0. Assum- 
ing that masses are smooth functions of the number of neutrons and pro- 
tons except for pairing effects, simple expressions have been derived for the 
gap parameter A based on binding energy differences between even and odd 
neighboring nuclei|0, [§. Detailed analyses^ and pairing adjustments^] have 
been based on these expressions. The simplest example is the well-known 
three-point mass formula: 

(-■\)N 

A^(N) = ^^[E(N+1)-2E(N) + E(N-1)] , (1) 
* Corresponding author : duguet@spht.saclay.cea.fr 
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where N is the number of nucleons (neutrons or protons). 

A study of the OES in light alkali-metal clusters and in light N = Z 
nuclei || has led to the conclusion that this phenomenon was not due to 
pairing correlations but rather to deformation effects (Jahn- Teller OES) [J7L 
[§[. That work motivated a study by Satula et al. || on the mean-field 
contribution to the OES in nuclei, especially coming from deformation. To 
isolate mean-field effects, the pairing force was set to zero and Hartree-Fock 
(HF) calculations were performed for light deformed nuclei. In this context, 
they observed an OES of the energy through A^ 3 \ itself oscillating: 



Aj&(iV) « if N is odd, 

~ 2 1S eVen • 

For even N, (ejt — et-i) is the gap around the Fermi level in the single- 
particle spectrum. It is zero for spherical nuclei (apart across sub-shells) and 
differs from zero for deformed nuclei because of the spread doubly-degenerate 
spectrum (see Fig. [l]). Deformation is thus found to be responsible for a direct 
contribution to the three point odd-even mass formula. 

With pairing correlations, A^ 3 '(odd) can a priori be a measure of pairing 
effects only, whereas A( 3 )(et>en) contains an additional contribution related 
to the splitting of the single-particle spectrum around the Fermi level. Such 
a scheme cannot account for the same oscillation of A^ 3 ^ in spherical nuclei 
because of the large degeneracy in spherical shells. 

On the other hand, recent calculations of spherical tin isotopes including 
pairing (HF+BCS calculations) JjlJ have led to the conclusion that the five- 
points formula A^ 5 ^ was in this case a better approximation of the pairing 
gap that A®(odd). 

The purpose of the present study is to analyze, for spherical and de- 
formed nuclei, the different contributions to odd-even mass differences in a 
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Figure 1: Schematic single-particle spectra and occupations for three succes- 
sive deformed nuclei. 

fully self-consistent mean-field picture including time-odd components of the 
force. Our aim is to give a coherent picture and to extract a quantity more 
directly related to pairing correlations. In addition, we want to investigate 
the connection between results obtained with and without the inclusion of 
pairing (cf. Eq. |2|), as well as the physical content of this connection. 

The present work is based on the conclusions of Ref. |TTJ] (hereafter ref- 
ereed as paper I), and is organized as follow. In section [| we introduce a 
separation of the microscopic binding energy which allows to separate dif- 
ferent types of contributions to odd-even mass formulas. The theoretical 
framework used to perform our mean-field calculations is detailed in sec- 
tion |3]. Results on spherical as well as deformed nuclei are presented in 
section f|. The evolution of odd-even mass differences as a function of pair- 
ing correlations intensity is studied in section [| using a schematic BCS model 
and self-consistent calculations. Finally, the analysis of the results and the 
conclusions are given in section ||. 



4 



2 Odd-even mass differences in self-consistent 
mean-field calculations 



To evaluate and understand the contributions to the odd-even mass differ- 
ences of nuclei in a fully self-consistent mean-field picture, two questions are 
addressed in what follows: 

(1) how to define a procedure to extract different contributions to the OES 
and to identify unambiguously their physical content? 

(2) is the analysis of the odd-even mass differences at the HF level of any 
help to understand what happens in presence of pairing correlations? 

2.1 Smooth contribution to mass formulas 

Several finite-difference mass formulas |3], f|, Q are used to evaluate the neu- 
tron or proton "pairing gaps". The aim is to extract the quickly varying 
part of the energy as a function of some parameters such as the number of 
neutrons or protons. The underlying assumption is that the microscopic en- 
ergy splits into a quickly varying part and a smooth one. In the description 
of nuclear structure, the rapidly varying component of the energy can be 
related to different phenomena such as shell closures, N = Z line, light mass 
nuclei, time-reversal symmetry breaking and reduction of pairing by block- 
ing in odd nuclei. The OES being related to the last two effects, appropriate 
mass regions must be chosen in order to avoid the first three ones. 

We define the smooth part of the energy as the one obtained when all 
nuclei are calculated as if they were even ones (no blocking and no breaking 
of time-reversal invariance in odd nuclei). Such an energy should not un- 
dergo odd-even irregularities. It will be referred to as ~Q HFBE for "Hartree- 
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Fock-Bogolyubov Even" and the associated wave-function will be denoted 
| ^,hfbe > Such a definition of the smooth part of the microscopic bind- 
ing energy has already been used in a work dealing with the OES in odd 



nuclei [10]. Then, the energy of an odd nucleus can be written as: 



E HFB (N) = E HFBE (N) + [E HFB (N)-E HFBE (N)} 



= E HFBE (N)+ E pol (N) + A (AT) , (3) 

where E pol (N) is the difference of binding energy due to the time- reversal 
symmetry breaking effect, and A(iV) is the positive contribution related to 
the fact that the odd nucleon is unpaired in the final HFB one quasi-particle 
(qp) state. It is denoted as the self-consistent pairing gap whereas E pol [N) + 



A(iV) is the self-consistent qp energy [jTT 



We shown in Ref. |Tl| that the state | \]/ HFBE > arises naturally as an 
intermediate step in the nucleon addition process. It defines an underlying 
even structure in an odd nucleus. 

2.2 Finite difference mass formulas 

Starting from the separation procedure defined by Eq. |, the smooth part of 
the binding energy can be expanded in a power serie around a given mass 
number N Q : 



oo i f)k rpHFBE 

E HFB (N) = ]T 



fc=0 k\ dN* 



(N - N ) k + E po \N) + A(N) (4) 



No 



Finite-difference formulas have been derived || f|, || to eliminate the 
successive derivatives of the smooth part of the energy. The three-point 
difference is written as: 
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^■hfb( n o) - Ahfbe(Nq) + Ap airing+pol (N ) . (5) 
Using Eq. |, we have: 

(-l) N o QlfiHFBE 



2 ON 2 
and 



(6) 

N 



&?Ln 9+P ai(No) = A(N ) + Ev° l (N ) ifiV isodd, 

= i{A(iV -l) + A(iV + l) 

+E p° 1 (JV - 1) + EP° l (N + l)} if N is even. 

Similar expressions are obtained for the fourth-order formula (five points 
difference) : 



( — 1 \ N o 

A<g> FB (N ) = -LJJ_ [E hfb {Nq+2) _ 4E hfb {No + 1) 

+ QE HFB (N ) - AE HFB (N -1) + E HFB (N -2)] (7) 

, (5) (AT \ I A ( 5 ) 



A HFB e( N o) + \airing+pol( N o) ■ 



Higher order formulas can be derived in the same way. 

2.3 Remarks on E HFBE and A{N) + E pol (N) 

If E HFBE represents a smooth part of the binding energy, its contribution to 
AW(N) should be of the same order of magnitude (at least in absolute value) 
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for odd and even neighbor nuclei. For instance, the contribution to A HFB (N) 
is approximately equal to {—\) N /2d 2 E HFBE /dN 2 , and since this second 
derivative is always positive and smooth, it gives an oscillating contribution 
of similar amplitude for odd and even N. 

Moreover, if such an hypothesis is valid, the contribution of E HFBE to 
A (n) (iV) will tend to zero, and A^ FB (N) to A^J iring+pol (N) , with increasing 
order n of the finite difference formula. The decrease of A^ BFBE with n will 
have to be checked in order to validate the above energy separation. 

With increasing n, A^ iring+ ol (N) is an average of A + E po1 over a larger 
number of nuclei. This assumes that A + E po1 varies slowly with the (odd) 
nucleon number. If so, A^ iring+pol (N) will be a relevant observable in the 
nucleus with iV nucleons. 

The choice of the order n of the formula used to extract a "pairing gap" 
A (AT) through the odd-even staggering should be a compromise as regards 
to the last two remarks. 



3 Theoretical framework 

In order to study the odd-even mass staggering in the framework of self- 
consistent calculations, one has to use the Hartree-Fock-Bogolyubov (HFB) 
approximation, since time-reversal symmetry is lost in odd nuclei because of 
the blocking of a single nucleon. This implies that time-odd components of 
the effective interaction must be included in the calculations. They are not as 
well determined []T^| as the even components since nothing constrains them 
specifically in the usual fitting procedure of effective forces [O . Despite these 
uncertainties, it is important to study their contributions to binding energies 



and odd-even mass differences which are known to be significant [|14| , p~5|j . 



The general formalism used here is detailed in Ref. |I6|. It is based on the 
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self-consistent Hartree-Fock Bogolyubov method, with an approximate parti- 
cle number projection by the Lipkin-Nogami prescription. In the particle-hole 
channel, we use a two-body force of the Skyrme type, SLy4, which has been 
adjusted to reproduce also the characteristics of the infinite neutron matter 



and, consequently, should have good isospin properties [13 . This force has 



been shown to describe satisfactorily nuclear properties for which it had not 



been adjusted such as super-deformed rotational bands [[17], and the struc- 
ture and decay of super-heavy elements Jl9[ . It should be mentioned that the 
time-odd components of the force are deeply involved in the description of 
rotational properties. The capacity of SLy4 to reproduce these observables 
is an advantage as regards to the previous discussion. In the T = 1 particle- 
particle channel, we use a surface-peaked delta force (Eq. H) adjusted on the 
low spin behavior of the moments of inertia of super- deformed bands in the 



A m 150 region | I7| . This pairing force has also been shown to work well in 



very different mass regions, up to the transfermium one, to describe ground 



state as well as rotational properties [20 . It is given by 



K = ^(l-P CT )5(ri-r 2 )(l-^), (8) 

* Pc 

where V T = — 1250 MeV.fm" 3 (r stands for proton or neutron), p(R) is the 
local matter density of the nucleus, P a is the spin exchange operator and p c 
= 0.16 fm -3 the nuclear saturation density. As V T is a contact interaction, 
we use a cut-off for the active pairing space which ranges from 5 MeV below 



to 5 MeV above the Fermi level [21] in the single-particle spectrum. 



4 Results 



9 



oA 3 


HFB 




HFB 



V 



•<> o Vo 



OA 


EXP 


• A 


EXP 



Figure 2: Left: calculated odd-even mass differences A^ B and A^ FB for 
the tin isotopic line from 100 Sn to 169 Sn. Right: known A^L and A^J. . 



Experimental data are taken from [22 



4.1 Spherical nuclei 

Seventy ground-states have been calculated along the tin isotopic chain, from 
100 Sn to 169 Sn. Each odd N nucleus has been calculated twice: first, as a HFB 
fully paired vacuum with an odd average number of neutrons (HFBE state) , 
then with the fully self-consistent HFBLN scheme. In this case, several one qp 
configurations are investigated to determine the configuration corresponding 
to the ground state. Due to the magic number of proton [Z = 50), all these 
nuclei are found to be spherical in our three-dimensional lattice calculation. 

Calculated and experimental p2| A^ 3 '(N) and A^(N) along this chain 
are given on Fig. [2| The staggering of A( 3 )(iV) can be seen on the two 
panels, the values for odd N being smaller than for even N. The amplitude 
of this staggering is smaller for neutron-rich nuclei. No staggering occurs for 
A^(N). 
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Figure 3: Calculated odd-even mass differences A#^ B£ and A^ BFBE for the 
tin isotopic line from 100 Sn to 169 Sn. 

Calculated A^ 3 ' 5 \N) are greater than experimental ones by several hun- 
dreds keV while the fit used for the pairing force has been shown to work well 
in reproducing rotational properties in several mass regions. This is partly 
due to numerical difficulties to converge very well several odd nuclei in the 
region 107 ~ 125 Sn which leads to an underestimation of their binding energy 
by two to three hundreds keV. Moreover, correlations beyond mean-field in 
odd nuclei should be larger than in even ones because of the large number 
of low-lying individual excited states in the first case. Treating explicitly the 
residual interaction through configuration mixing in odd and even nuclei is 
expected to lower the OES by approximately three hundreds keV |23| . 



On Fig. H are plotted the contributions A BFBE and A BFBE to A BFB and 
^hfb respectively. Apart for the magic number N = 82, one has |A^ BB | ^> 
\A BFBE \ = 0. This justifies the identification of E HFBE as the smooth part 
of the energy, and indicates that higher order formulas are not needed. 
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Figure 4: Left: A HFBE is compared to A HFB - A HFB along the tin isotopic 
line. Right: comparison of A HFBE with experiment (see text). 

Assuming that A® iring+pol is equal to A { ^ iring+pol since A(N) + £ poi (iV) 
is constant over a few odd nuclei (see section |2.1| ), one can write: 



A§ FB {N) - A® FB (N) » Ag } F ^(iV). (9) 

The two sides of Eq. |9| are compared on the left handside of Fig. [|. The 
accuracy of the agreement is impressive along the whole isotopic line (except 
as expected for the magic number N = 82). It confirms that the contribution 
of the smooth part of the energy to A^ 3 ) is equal in absolute value for even and 
odd N (see section |2?l| ), the sign (— 1)^ being responsible for the oscillating 
pattern of this contribution. The A^ 3 ) staggering corresponds to an oscillation 
around A^ 5 ) due to A ( F[ ] FBE . One can therefore conclude that A^ 5 ^ is a measure 
of the rapidly varying part of the energy A^ iri + AN) « A (AT) + E pol (N). 
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Figure 5: Left: calculated odd-even mass differences A H f FB and A^ BFB for 
the tin isotopic line with a reduced neutron pairing strength V n = — 1000 
MeV.fm -3 . Right: A^ FBE is compared to - A^ FXP for the same V n . 



On the contrary, A^(odd) contains smooth contributions, in particular the 
full asymmetry energy contribution to the OES. 

The right panel of Fig. |] provides a comparison between A BFBE and 
Afil p - Ajp x p . The agreement is very good. Although the absolute values 
of Ag^p are overestimated by a few hundred keV, the staggering of 
around A^j. is well reproduced by the contribution coming from the smooth 
part of the energy A HFBE . This shows the strong decoupling between the 
two contributions to Eq. ||. 

To confirm this decoupling, we have performed the same calculation with 
a decrease of the neutron pairing intensity V n by 20 % to —1000 MeV.fm" 3 . 
The left panel of Fig. |5| shows the response of the absolute OES A BFB . 
A few points are missing on the figure due to problems of convergence for 
some nuclei with a reduced pairing strength. The odd-even mass differences 
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A-hfb an d ^hfb which include the contribution coming from the pairing 
gap, decrease also by approximately 20 %. On the other hand, the right 
panel of Fig. [5] shows that the oscillation of A^ FB due to the smooth part of 
the energy is not modified. Thus, the separation of Eq. |3| divides the energy 
in a part extracted through responding directly to the pairing intensity 
and in another one coming from E HFBE (almost) insensitive to it. 

In a previous study of 254 No PU|, we have shown that the dynamical 
moment of inertia of its ground state rotational band depends more strongly 
on the radial dependence of the pairing force than on its intensity for a given 
radial form factor. On the contrary, the OES around 254 No has been found to 
vary proportionally to a change of the intensity of the pairing. From this and 
from the present results, we can conclude that the part of the pairing energy 
contained in the even energy E HFBE is probed by observables involving the 
nucleus as a whole, such as rotation, and is more sensitive to the analytical 
structure of the force. On the other hand, the part defined by A(iV) is related 
to a specific blocked orbit and is very sensitive to the intensity of the pairing 
force. This provides with two different observables to adjust the strength 
and the radial dependence of the pairing force. 



4.2 Deformed nuclei 

Let us now extend our analysis to deformed nuclei. Forty nine Cerium iso- 
topes have been calculated, from 118 Ce to 166 Ce. The deformation parameter, 
/?2, is given as a function of A on Fig. ^| for HFB and HFBE calculations. 
The ground-state quadrupole deformation undergoes large variations from a 
region of strong prolate deformation around 118 Ce (/?2 ~ 0.37) to a prolate 
deformation region (/5 2 ~ 0.31 around 160 Ce) through the spherical 140 Ce nu- 
cleus. Fig. H illustrates the fact that $2 FBE reproduces the mean evolution 
of /3f FB with A. 



14 



P 2 u 



1 / 



• HFB 
« HFBE 



120 130 140 150 160 170 



Figure 6: Quadrupole axial deformation as a function of the mass number 
along the cerium isotopic chain for HFB and HFBE calculations. 

It has been shown in paper I that, apart for the time-reversal symmetry 
breaking effect, the HFBE and HFB calculations lead to nearly identical 
mean-fields. The influence of the deformation on the binding energy of odd 
nuclei is contained in E HFBE and the Jahn Teller contribution to the odd-even 
mass differences ^hfb wm thus be extracted through A BFBE . The energy 
separation [3] should allow to isolate the self-consistent qp energy through odd- 
even mass formulas in the same way as for spherical nuclei. Let us illustrate 
these statements along the cerium isotopic line. 

On the left (right) panel of Fig. [F] are plotted A BFB ^ and Ajj FB , E y 
The comparison between A BFB — A BFB and A BFBE is presented on Fig. |[ 

(3) 

Similar results and agreements as for tin isotopes are obtained: \A HFBE \ is 
much larger than \A BFBE \ which is close to zero and A H i FB (N) — A BFB (N) 
is approximately equal to A BFBE (N). The contribution A(N) + E pol (N) 
to the odd HFB energy is extracted through A^ for deformed nuclei as for 
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Figure 7: Left: calculated odd-even mass differences &hfb and A HFB for the 
cerium isotopic line from 118 Ce to 166 Ce. Right: calculated odd-even mass 
differences £±hfbe an< ^ ^hfbe- 

spherical ones. The effect of deformation on the OES (Jahn Teller OES) 
is included in the contribution coming from E HFBE which reproduces the 
oscillation of A^ 3 ) around A^ 5 ) (Fig. |8|). This contribution is identical for 
odd and even neighboring nuclei, with an opposite sign contrary to what was 
found in Ref. ||. We will come back to this point later. It is unfortunately 
not possible to make a significant comparison with the experimental data 
for these deformed nuclei because the experimental error bars are much too 
large. 



5 Transition to zero pairing 

Two types of calculations have been performed to study the OES in the zero 
pairing limit. We have first used a schematic non self-consistent BCS model 
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Figure 8: Comparison between calculated A HFBE (N) (filled circles) and 
A^ FB (N) — A BFB (N) (empty diamonds) along the cerium isotopic line. 

in order to get qualitative informations. Then, fully self-consistent HFB 
calculations have been performed in order to take into account rearrangement 
effects and obtain quantitative results. 



5.1 Schematic model presentation 

The schematic model consists in a non self-consistent BCS scheme on top 
of a fixed single-particle spectrum {e^} for one kind of particles only. No 
two-body force is included: the orbital dependent gap is given as an initial 
parameter and is parametrized as: 



A fc = A Exp 



3A 



(10) 



A being the input. Let us mention that A is the gap at the fermi energy. 
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The fixed spectrum can be either a priori constructed or taken from the 
self-consistent calculation of an even nucleus as a typical spectrum in a nar- 
row region around that nucleus. In what follows, the subscripts HFB(E) are 
changed into BCS(E). 

The fully paired part of the energy as defined by the first term of Eq. |3] 
is that of a fully paired BCS vacuum: 

e bcse {n) = J>2 e * £ Ag , (11) 

k 2 k>o J(e k - X N ) + A\ 



X N being fixed by the additional condition < N > = N (odd or even). For 
N even, it corresponds to the exact BCS energy. For odd N the BCS energy 
is: 



E BCS (N) = E BCSE {N) + Mm {Ef} , (12) 

where the Ef(N) = J (e^ — A^) 2 + A| are the quasi-particle energies eval- 
uated in the odd vacuum JET]]. Once these energies are given, the energy 
differences A^ cs , A^ CSE , A^ cs and A^ CSE can be computed. 



5.2 Results on cerium isotopes 

We have performed two different applications of the schematic model. First, 
we have used an equidistant two-fold degenerate spectrum simulating a de- 
formed nucleus with a single-particle level spacing 8e = 400 keV. The calcu- 
lation has been performed for six different values of A, from zero to a typical 
value of 1200 keV. In the second case, we have used a realistic spectrum 
of 152 Ce (self-consistent HF spectrum obtained with the SLy4 interaction). 
Fig. |9] displays the evolution of Ag CS , A^ cs with increasing gap A from 148 Ce 
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Figure 9: A^ s , A^ cs and A as a function of A. From top to bottom, the 
pairing gap increases from to 1200 keV. Left column: calculation with an 
equidistant doubly degenerate spectrum. Right column: calculation using 
the self-consistent neutron HF spectrum of 152 Ce. In addition to A, the 
lowest qp in odd nuclei is shown. Results are displayed between 148 Ce and 
156 Ce. 
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to 156 Ce for both spectra. The left column is for the equidistant spectrum 
while the right one is for the realistic spectrum. 

For A = 0, we recover qualitatively the results obtained in HF calculations 
without time- reversal symmetry breaking (cf. Eq. |2|). Namely, A B ^ CS (N) 

(5) 

oscillates between for odd N and a non zero value for even N; A B ^ S is 
different from zero for all values of N. It follows that in this case, A^g^odd) 
extracts the gap. 

When A is increased, this is no longer true. Odd-even mass differences 

(5) 

are shifted to higher values in such a way that A B ^ S tends to extract the 
pairing gap while A B c S oscillates around this value. From the equidistant 
spectrum calculation, one can see that this asymptotic situation is achieved 
for a ratio A/5e ~ 0.5; i.e as soon as A reaches 200 keV. In real nuclei, apart 
for near closed-shell nuclei, typical A/5e values are greater than 0.5. 

The right hand side of Fig. shows that similar qualitative results are 
obtained using a realistic spectrum for deformed cerium isotopes, even though 
the structure of the realistic spectrum modifies the artificial regularities of 
the former case. The discrepancy between A and Min {Ef} observed for 
some nuclei is due to the fact that the qp in the odd fully paired vacuum 
is not always such that u 2 — v 2 is exactly JTT] . From a quantitative point 

(5 s ) 

of view, A B lj S extracts precisely the quasi-particle energy as soon as A at 
the Fermi energy reaches 60 % of the realistic value obtained in an HFB 
calculation of these isotopes. 

This can be understood from Fig. [iy which gives the E BCSE contribution 
to the different odd-even mass formulas for identical values of the gap, using 
both spectra. The oscillation of A BCS (N) around A BCS (N) is also plotted. 
For A = 0, A^cse is equal to A^ B ^ S as in this case E BCS and E BCSE are 
the same for odd nuclei as well (Min {E^} = in the odd vacuum for a 
vanishing pairing). Through the even part of the energy E B , we have 
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Figure 10: Energy differences ^-bcse^ ^bcse an d ^bcs " ^bcs- Left and 
right panels: same as for Fig. |9|. The results are displayed between 148 Ce and 
156 Ce. 
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isolated the quantity responsible for the Jahn Teller OES of Eq. ^ in absence 
of pairing. However, as the pairing increases this part of the energy becomes 
smoother with N, in such a way that A ( § > CSE goes to zero while A B ^ CSE 
oscillates regularly around it. When pairing increases, the deformation effect 
on the OES is equally redistributed to even and odd A^(N) in such a way 
that it produces the staggering of around Ag CS . 

These schematic calculations help understanding the link between the 
apparently contradictory results obtained in HF calculations where is 

(5) 

not zero and the BCS ones with the pairing turned on where A BCS equals 
Ap a i r i n g( + p i). The contributions to A^ 5 ^ coming from the two parts of the 
energy evolve in opposite ways. The even contribution decreases from a non 
zero value because of the deformation effect described in Ref. [|], to zero 
value with increasing pairing. At the same time, the blocking contribution 
Apairing(+poi) increases with pairing as expected. This illustrates why and 
how A^ 5 ) extracts A pairing ( +po ^ only, for a realistic pairing strength. 

Fig. |ll] gives a graphic representation of E HFB and E HFBE as a function 
of TV for both the zero pairing and the realistic pairing cases. It is drawn 
for an underlying doubly degenerate single-particle spectrum typical of a 
deformed nucleus. The left panel gives the zero pairing case where one can 
see that the even part of the energy is not a smooth function of iV and that 
the exact odd energy differs from E HFBE only if polarisation is included. The 
Kramers degeneracy is responsible for the linear behavior of E HFBE between 
two successive even nuclei while the single-particle level spacing is responsible 
for the change of the corresponding slope and induces the asymmetric Jahn 
Teller OES of Eq. ^. In the realistic pairing case displayed on the right 
panel, the behavior of E HFBE becomes smooth with the nucleon number and 
no asymmetry between odd and even N remains. As a consequence, one 
can graphically see that the Jahn Teller OES cannot be transposed to the 
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Figure 11: Binding energy as a function of N for the even part (squares joined 
by dashed line) and for the full odd states (circle). Left panel: no pairing. 
Right panel: realistic pairing case. 

realistic pairing case as it is in the absence of pairing. 
5.3 Results on tin isotopes 

For spherical nuclei, the same kind of A^ 3 ) staggering as for deformed nu- 
clei is observed experimentally and found theoretically in HFB calculations, 
while such staggering does not occur at the HF level because of the strong 
degeneracy of the spherical shells. 

We now apply our model with the HF spectrum of 122 Sn. Fig. [12| displays 
the same quantities as Fig. ^] between 118 Sn and 126 Sn, with a gap varying 
between and 1900 keV. This latter value corresponds to the theoretical 
gap at the Fermi energy in the HFB calculations. As expected, no Abcs 
staggering is seen for A equals apart for the transition between 2d3/2 and 
lhll/2 spherical shells which occurs at N equals 120. 



23 



> 

0> 



600 
500 
400 
300 
200 
100 

800 
700 
600 
500 
400 

1500 
1400 
1300 
1200 
1100 

2100 
2000 
1900 
1800 
1700 




I T I i I i I i I 




I I I I I I I I I 




_l i I I I I I i L_ 



A BCS 
*■ ■ A 5 BCS 

^«MinE 

qp 

A^A A = OkeV 



A= 400 keV 



A— A A= 1200 keV 



A— A A= 1900 keV 



118 120 122 124 126 

A 

Figure 12: Same as Fig.|]for tin isotopes. The calculation is done using 122 Sn 
HF spectrum. Results are displayed between 118 Sn and 126 Sn. 
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As A increases (from top to bottom on Fig. [12|), two modifications on odd- 
even differences occur simultaneously, namely the appearance of an odd-even 
staggering for and the extraction of Min {Ef} through A^ cs . As for 

deformed nuclei, the oscillating behavior of A^ s is directly related to the 
contribution from the even part of the energy. This transition from a situ- 
ation where no ^bcs staggering exists for zero pairing to a situation where 
a clear staggering develops, shows that the calculations with and without 
pairing are not related in a simple way. 



5.4 Single-particle level spacing and odd-even mass 
formulas 



In Ref. 0, |24]], it has been suggested that the A^ oscillation as a function of 
N could be used as a measure of the splitting of the single-particle spectrum 
around the Fermi level for even deformed nuclei (see Eq. ^). To study the 
validity of this statement, we compare on Fig. [13] A^ cs (2j) - A^ cs (2j + 1) 
and {ck+\ - efc)/2 for zero and a realistic value of the pairing strength. The 
variation of the chemical potential as a function of A is also plotted. In 
the left column are presented the results of the calculations for the cerium 
isotopes with an equidistant spectrum and with a realistic spectrum in the 
middle. In the right part of the figure are shown the results for tin isotopes 
using a realistic spectrum. 

For A = and 1200 keV, the splitting of the single-particle energies 
in the cerium spectrum is exactly reproduced by the staggering of A® if 
the equidistant spectrum is used. This result suggests that one can extract 
informations about the neutron (proton) single-particle spectrum along an 
isotopic (isotonic) line through odd-even mass differences whether pairing is 
present or not. On the other hand, using a realistic spectrum this conclusion 
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is valid only in the limit of vanishing pairing. For a realistic value of A, 
the difference A ( B ) cs (2j) - A^ cs (2j + 1) is no longer a measure of (e& + i - 
efc)/2. This is actually the case as soon as A reaches 60% of a realistic 
value. This is further confirmed by the calculation on tin isotopes where 
A BCS (2j) - A BCS (2j + 1) is non zero for realistic pairing strengths whereas 
the corresponding spherical single-particle energies are highly degenerate (see 
the right column). 

This result can be understood as a consequence of the very different way a 
nucleon is added whether pairing correlations are present or not. When start- 
ing from an HF state, the pairing is increased, the amount of binding energy 
associated with the addition of a nucleon in the even structure is less and less 
related to a specific single-particle energy. Rather, the nucleon is spread out 
on the levels around the Fermi level because of pairs scattering [p| . Con- 
sequently, the memory of the underlying single-particle spectrum is washed 
out. Besides, this is the reason why the HFBE energy becomes smoother as 
a function of A with increasing pairing (cf. Fig. |Tl|). 



The bottom row of Fig. [13] illustrates the previous statement by showing 
the chemical potential in the even state as a function of the mass number. 
Let us consider the two cases which make use of a realistic spectrum (middle 
and right panels of the bottom row). For A = 0, A is sensitive to the orbits, 
while for A rea u st i C it behaves more smoothly as a function of A and does not 
reflect the structure of the spectrum anymore. Thus, since one has 

Ag» s (2 J ')-AgL(2j + l)*^-^ (13) 



this observable is no longer directly related to the single-particle level spacing 
around the Fermi level for a realistic pairing strength. The strong influence 
of the single-particle levels structure is lost for a value of the gap smaller 
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Figure 13: Upper rows: A^ 3 ^(2j) - A^ 3 ^(2j + 1) compared to the splitting 
around the Fermi level in the even nucleus (e^+i - efc)/2 for two extreme 
values of the gap (top, A = 0; middle, A rea i istic ). The 1.5 MeV splitting for 
120 Sn is out of scale. Lower row: chemical potential in MeV as a function of 
A for the two values of A. 
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than typical splittings in the spectrum (Ae goes from 200 keV to 1 MeV in 
the studied cerium region). 

In the calculation based on an equidistant spectrum, the left panel of the 
bottom row illustrates why in this case one can still extract informations 
about the single-particle level spacings through odd-even mass differences 
for a realistic pairing intensity. The evolution of A with A does not depend 
on A. Indeed, even if the nucleon is spread over the Fermi sea, the average 
effect of the pair scattering process cancels out because of the symmetry of 
this spectrum. As a result, the energy added by the extra nucleons remains 
equal to the single-particle energy of the orbit on which the nucleon is put in 
absence of pairing. This points out the inherent limits of schematic models 
used with very simplified single-particle spectra. 

The above result could not have been worked out for spherical nuclei with 



models limited to a single j shell. Indeed, from the bottom row of Fig. [13|, one 
sees that the effect involves several spherical shells for the pair scattering. For 
instance, the pair scattering effect is efficient enough for A equal to A rea i istic 
to loose the information about the 1.5 MeV splitting between the 2d3/2 and 
lhll/2 spherical shells. 

5.5 Self-consistent calculations 

Let us present the same analysis for fully self-consistent HFB calculations 
of cerium isotopes. Pairing correlations are gradually turned on through 
the increase of the neutron pairing force intensity V n up to the realistic case 
presented in section P~2| . Fig. |TJ] displays the same quantities as Fig. [|and [10 



Instead of the perturbative BCS quasi-particle energy Min {Ef} , the energy 
difference E HFB - E HFBE is given. This quantity is the self-consistent version 



of the created quasi-particle energy in odd nuclei [11]. The upper left panel 



shows that it is non-zero in the zero pairing case since it already contains the 
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Figure 14: Left (right) column: A^ FB ^ (circles with solid line), ^hfb{e) 
(squares with dashed line), self-consistent qp energy E HFB - E HFBE (triangles 
with dotted line) and A BFB — A BFS (stars with dotted line). From top to 
bottom the intensity of the neutron pairing force increases from to 1250 
MeV.fm -3 . Results are displayed between 148 Ce and 156 Ce. 



time-reversal symmetry breaking effect which shifts up all the A^ BFB . 

First, let us concentrate on the odd-even differences of Yi HFBE for V n = 
(right upper panel). For this particular case, we will use HF and HFE sub- 
scripts instead of HFB and HFBE. There are neither polarisation nor pair- 
ing effects and one can focus on the strong influence on the results of the 
self-consistency in the mean-field treatment. Contrary to the schematic re- 
suits, A HFE [odd) can be significantly different from zero in this case (several 
hundreds keV for example between 145 Ce and 149 Ce). Self-consistency signif- 
icantly modifies the picture as compared to the independent particle scheme, 
especially in regions of varying deformation || p5] . Since A HFE (odd) is not 
0, it is thus difficult to argue that it extracts A pa i r i ng . The only possible 
statement is that in the zero pairing case an OES is seen with an oscillating 
Ajjp together with a non-zero A^ F . 

The energy difference A HF (odd) is closer to E^ - E HFE than A^p. How- 

(5) 

ever, as the pairing intensity increases A HFB extracts the energy difference 
•^hfb _ j^hfbe ^ -y^i^ j s nothing but the staggering of the energy associated 

(3) 

with blocking, while A HFB oscillates around this value. This is the case as 
soon as the pairing intensity reaches about 72 % of the realistic value (V n = 
— 900 MeV.fm -3 ). This statement is valid even when self-consistency effects 
are large in this region of varying deformation. The results are presented 
only for a small part of the cerium isotopic line, the same conclusions holds 
for the whole line. 

From the right column giving the even contributions for increasing strength 
V n , one can see that the analysis of the schematic calculations remain true. 
Namely, as also discussed in section [17^, the oscillations of A^ FB around 

(5) 

hfb are reproduced quantitatively by the contribution from the even, 
smooth, part of the energy for realistic V n . 

,(3) Q ^ A ( 3 ) 



On Fig. 15, the A HFB and A HFBE staggerings are compared with the 
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Figure 15: Energy differences A>£> FB (2j) - A% FB (2j + 1), A% FBE (2j) - 
A HFBE (2j + 1) and single-particle splitting (e^+i — ek)/2 around the Fermi 
energy for self-consistent calculations of cerium isotopes. 

single-particle level spacing around the Fermi energy in even nuclei. The left 
panel displays A% FB (2j) - A® FB (2j + 1), A® FBE (2j) - A% FBE (2 3 + 1) and 
(efc+i — e&)/2 for V n = 0. The staggerings of A HFB and A HFBE coincides. 
It means that this staggering is entirely due to the even contribution to the 
EOS. Moreover, the A (3) stag gering roughly extracts the splitting in the HF 
spectrum except in the region of varying deformation ( 124 Ce to 152 Ce) where 
the rearrangement due to self-consistency from one nucleus to the next is 
large. 

The right panel displays the same quantities in the case of realistic pair- 
ing intensity. With pairing included, e£"^° is the eigenenergy of the HF field 
deduced from e£ n ~° by continuity. Again, the A BFB staggering is well repro- 
duced by A BFBE (2j) — A BFBE (2j + 1) along the whole line, whatever the 
magnitude of the self-consistency effects is. However, the information about 
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the HF eigen-energies is lost in this case. Indeed, the addition of a nucleon 
is no longer related to a single orbit when pairing is included fil l. These 
conclusions are the same as in the case of the schematic model. 



6 Analysis and conclusions 

We have proposed an analysis for the odd-even mass staggering based on the 
definition of a "virtual" odd nucleus (HFBE state) having the structure of 
an even one as the underlying structure of the "real" odd nucleus |T0|, [TI| . 

For realistic pairing intensities, it has been shown that the A^(N) mass 
formula extracts precisely the self-consistent HFB quasi-particle energy for 
spherical as well as for deformed nuclei. The self-consistent HFB quasi- 
particle energy corresponds to the blocking of the odd nucleon on top of the 
fully paired odd reference vacuum and contains both the pairing gap, A(N), 
and the time- reversal symmetry breaking effect, W° l . 



Similar results have already been reported for spherical nuclei in Ref. [[10 
where the extraction through A^ 5 ^ of the pure blocking contribution to the 
ground-state energy of odd nuclei was pointed out. However, this work was 
done in the HF + BCS framework and without breaking time-reversal sym- 
metry and A^ 5 ) extracts then only the self-consistent pairing gap (W° l = 0). 
Such an approximation limits the pertinence of the comparison with exper- 
imental data. Our present study does incorporate this physical effect and 
extends that earlier work to realistic cases and to deformed nuclei. 

Satula et al. [[| made a similar study for deformed nuclei in the HF 
approximation, as a reference to identify pairing contributions to the OES. 
They proposed to use A^\odd) as a measure of the pairing gap and the 
difference A^ 3 \even) - A^ 3 \odd) as the Jahn Teller contribution to the OES 
(contribution from deformation). They have also neglected time-reversal 



32 



symmetry breaking effects. 

In this context of time-reversal invariance, our extended analysis of the 
OES as a function of pairing correlations intensity within the frame of a 
schematic BCS model has allowed to sort out the contradictory former propo- 
sitions. The conclusions of Satula et at, based on Eq. 0, have been shown to 
hold only for very weak pairing. For a realistic pairing intensity, A^(odd) is 
no longer a measure of the gap alone since it contains an additional contri- 
bution coming from the even part of the energy E HFBE as defined in Eq. |3] 
(see section [5]). On the other hand, A^ 5 ^ extracts the pairing gap in this case. 
Self-consistent HFB calculations have confirmed these conclusions from a 
quantitative point of view. 

Let us go one step further by introducing, in a self-consistent mean-field 
picture, the time-reversal symmetry breaking effect on binding energy. This 
effect is formally related to the physical blocking process in odd nuclei as 
extensively discussed in Ref. |TT]. It follows that it is deeply associated 
with the self-consistent pairing gap in such a way that these two energetic 
quantities cannot be separated through odd-even mass differences. They are 
both contained into A( n \N) at all orders in n when using experimental data. 
Consequently, one has to include this effect in realistic calculations in order 
to compare directly theoretical and experimental odd-even mass differences. 

Finally, we have identified in the present paper the physical content of 
A^ 3 ) and A^ 5 ^ in fully self-consistent mean-field calculations including realistic 
pairing: 



/i\iV q2 TpHFBE 

^hfb(N) « A(N) + E pol (N) + ^- 



(14) 



2 dN 2 

^hfb(N) « A(N) + E pol (N). (15) 



N 

.(5) 
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where in the picture of Ref. [|TT|, 1/2 d 2 E / dN 2 is related to the nucleoli 
addition process and contains the full assymetry energy contribution to the 
OES whereas A(N) + E pol (N) is related to the blocking of this nucleon. 

Comparing their results with those obtained by Satula et al. ||, Bender 
et al. [|1(| argue that the Jahn Teller effect (called "mean-field effect" since it 
is related to the structure of the single-particle spectrum) is not connected to 
the oscillation of A^ 3 ^ found for spherical nuclei (E HFBE contribution). We 
have demonstrated that the A^ 3 ) staggering was always related to the E HFBE 
contribution (typically ±50/150 keV). This energy is related in some way to 
the single-particle structure of a given nucleus, but our extended schematic 
and fully self-consistent calculations have shown that the experimental A® 
staggering cannot be identified with single-particle level spacing at the Fermi 
surface as suggested in Ref. ||, apart for nuclei immediately near magic 
ones. 

Our results are based on calculations done in the A = 100-170 mass region. 
They should also be valid for lighter nuclei. Indeed, the regime (independent 
particle or correlated system) in which the system stands depends on a typical 
ratio A/5e. It has been shown in a schematic BCS model that the correlated 
regime is achieved for a ratio A/5e ~ 0.5 whereas for realistic calculations in 
the A = 100-170 mass region, it is achieved for a value of the pairing gap a 
few times smaller than the level spacings near the Fermi energy. These two 
arguments are in favor of the correlated limit for nuclei in the mass region A 
= 30-100 where Ap/Sep is typically between 0.5 and 1 for mid-shell nuclei. In 
order to check this statement, we have performed an exploratory calculation 
for Mg isotopes between 24 Mg and 28 Mg using our schematic BCS model. 
The single-particle spectrum and gap value at the Fermi energy taken from 
a self-consistent HFB calculation of 24 Mg were used. The results support the 
extrapolation of our results to lighter masses. It also shows that the E HFBE 
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contribution to A^ 3 ^ increases in average with decreasing mass number as the 
mean single-particle level spacing increases at the same time. It qualitatively 
explains the well known increase of the A^ 3 ^ staggering around A^ 5 ) in light 
nuclei (see Fig. 3 of Ref. || for instance). 

The only limitation of the above conclusions concerns nuclei with neutron 
(proton) number one or two units away from magic numbers. These nuclei 
belong to the intermediate regime where A^(odd) is of the same quality or 
better than A^ 5 ^ to extract informations about the blocking effect (see upper 
panels of Fig. ^j. However, the limitation concerns a very limited number of 
nuclei which, in any case, should not be used for a study intended to adjust 
the pairing force. It has also been shown that the pairing force should be 
fitted on global observables such as rotation, and orbit-related observables as 
the OES in order to adjust at the same time the parts of the pairing energy 
contained in E HFBE and in A(7V). 

The above analysis is directly related to the nucleon addition process 
which is significantly modified by the inclusion of pairing correlations in the 



nuclear wave-function | IT| . Besides, in a recent lecture [25J where Satula 



and coworkers results on the OES were reported, Flocard suggested that 
it is somewhat surprising that the prescription of Eq. [2] derived using an 
independent particle picture remains correct for strongly correlated systems 
as nuclei. The present work has shown that this doubt was justified since 
pairing in such systems is strong enough in general to modify the picture 
by washing out the decisive influence of single-particle energies on odd-even 
effects. 

Once we have identified the physical content of the and odd- 
even mass formulas, their respective advantages and drawbacks as suited 
quantities for a detailed study or fit of a pairing force remain to be analyzed. 
Eq. [IR shows that A^ 5 ^ contains one quantity in addition to the pairing 
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Figure 16: Comparison between A HFB (odd) and the self-consistent qp en- 
ergy A + E po1 along the cerium isotopic chain. The comparison between 
Aff FB (odd) and an approximation of the self-consistent pairing gap A is also 
given. 



36 



gap. Actually, A^ 5 '(odd) for example contains a weighted average of A + E po1 
over three odd nuclei. As shown on Fig. |16| this is responsible for a slight 
deterioration of the validity of identity when A + E poX changes suddenly 
around one nucleus (see 163 Ce for instance). On the other hand, A^\odd) 
contains A + E po1 from the studied nucleus only which is an advantage over 
A( 5 ). 

Eq. [14] shows that A® contains an extra contribution coming from the 
smooth part of the energy E HFBE . In section |], this extra contribution 
has been shown to be of the order of ±50 to ±100 keV in spherical tin 
isotopes and of the order of ±100 to ±150 keV in cerium deformed nuclei, 
namely it contributes for about 8 to 12 %. Then, the time-reversal symmetry 
effect has been theoretically extracted in Ref. |l i] through a perturbative 
calculation (labeled Approx(E poi ) in the present work) for the cerium isotopes 
and appeared to be of the order of ±100 to ±150 keV. 

It follows that the two last previous contributions roughly cancel out in 
A^(odd) and that the relative weight of A(N) is larger in A^(odd) than 
in A^(odd). However, the details of this cancelation is not under control 
since A^ FBE and above all E pol (N), are not well known. In particular, 
the time-reversal symmetry breaking process deserves more studies since the 
results are force and model dependent [?, [H], 2J, 27]. In order to exemplify 
the situation, Fig. [16] also gives A^(odd) and an approximation of the self- 
consistent pairing gap, [A ± E po1 - Appiox(E po1 )]. Results are not shown 



when the hypothesis of the perturbative calculation are not fulfilled [|Tl 
One can see that A^ BFB (odd) is often closer to the self-consistent pairing gap 
than A^ BFB (odd) which means that the cancelation between the two different 
effects E pd and A^ FBE is quite effective in the present case. 

Finally, one should propose A^^odd) as the better suited quantity for a 
detailed study of the pairing gap or for the fit of a pairing force through the 
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adjustment of theoretical and experimental odd-even mass differences. We 
would like to stress the fact that this conclusion is not a validation of the anal- 
ysis performed in || as the way to reach it has been very different and needed 
the inclusion of time-reversal symmetry breaking in order to point out the 
a priori unexpected cancelation between E po1 and A^ FBE (odd). Moreover, 
this conclusion still depends on more extensive analysis of the time-reversal 
symmetry breaking contribution in different mass regions to be done in order 
to study the presently found cancelation effect. As an example, we may in- 
terpret this effect to be responsible for the overall agreement found between 
and several average d[] theoretical HFB pairing gaps in an extensive re- 
analysis of the commonly accepted formula A = 12 A^ 1 ^ 2 MeV for the pairing 



function of the mass number [E8I 



It is important to stress that our purpose takes into account only one kind 
of pairing correlations, i.e. proton-proton and neutron-neutron pairing. The 
questions related to proton-neutron cooper pairs around N = Z nuclei need 

as 



of course an extension of our approach. Satula and Wyss p9f, Vogel [30 



well Terasaki et al. |3T| have investigated these questions. Their conclusions 
correlated to an extension of our work could deliver a good indicator to fix 
the theoretical intensity of this neutron-proton pairing. 
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